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^ Abstract 

o\ 

,_i We present three sets of results for the stationary distribution of a two-dimensional 

semimartingale reflecting Brownian motion (SRBM) that lives in the nonnegative 
Oh quadrant. The SRBM data can equivalently be specified by three geometric objects, 

an ellipse and two lines, in the two-dimensional Euclidean space. First, we revisit 
the variational problem (VP) associated with the SRBM. Building on Avram, Dai 
^ and Hasenbein (2001), we show that the value of the VP at a point in the quadrant 

is equal to the optimal value of a linear function over a convex domain. Depending 
on the location of the point, the convex domain is either P^^) or P^^) pCi) nP^^) 
'i> where each T>^'^\ i = 1,2, can easily be described by the three geometric objects. 

Our results provide a geometric interpretation for the value function of the VP and 
allow one to see geometrically when one edge of the quadrant has influence on the 
optimal path traveling from the origin to a destination point. Second, we provide 
a geometric condition that characterizes the existence of a product form stationary 
distribution. Third, we establish exact tail asymptotics of two boundary measures 
,__( that are associated with the stationary distribution; a key step in our proof is to 

sharpen two asymptotic inversion lemmas in Dai and Miyazawa (2011) that allow 
one to infer the exact tail asymptotic of a boundary measure from the singularity 
of its moment generating function. 



1 Introduction 



Multidimensional semimartingale reflecting Brownian motions (SRBMs) have been ex- 
tensively studied in literature. They serve as diffusion models of multiclass queueing 
networks. See [1 1] for a survey on SRBMs and [ ] on diffusion approximations. In this 
paper, we focus on a two-dimensional SRBM Z = {Z{t);t > 0} that lives on the state 
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space M^. We adopt the definition of an SRBM used in Dai and Miyazawa [ ] and follow 
the notational system there. In particular, the data for the SRBM consists of a covariance 
matrix E, a drift vector /x, and a reflection matrix R. They are written as 

In short, Z is said to be a (E, /i, i?)-SRBM on M^. The covariance matrix E and drift 
vector II uniquely determine the following ellipse passing through the origin in M^: 

where {x,y) denotes the usual inner product for vectors x,y G (see Section 2 for 
more discussion on this inner product). The first column R^ of the reflection matrix 
R determines a line, called line 1 in this paper, that passes through the origin and is 
orthogonal to direction R^. In addition to the origin, the line intersects the ellipse at 
point 9^'^''^^ 7^ 0. Similarly, the second column R^ of R defines line 2 with intersection 
^(i-r) ^ on the ellipse. These three geometric objects, the ellipse and two lines, uniquely 
determine the SRBM data (S, /i, R). See Figure 1 for an example of these three geometric 
objects and Section 3 for a precise definition of these objects. All results in this paper 
will have a geometric interpretation in terms of these three objects. 

Throughout this paper, we assume conditions (2.1) and (2.2) in Section 2 on the SRBM 
data. These two conditions are necessary and sufficient for a two-dimensional SRBM to 
have a unique stationary distribution. Section 3.1 will give a geometric interpretation of 
these two conditions. In this paper, we present three sets of results for a two-dimensional 
SRBM. We first give a geometric interpretation for the value function of the variational 
problem that is associated with the SRBM. The value function has been conjectured 
to be the large deviations rate function of the stationary distribution. See inequalities 
(2.6), (2.7), and the text immediately following these inequalities for discussion on large 
deviations rate functions. We next give a geometric condition for the SRBM to have a 
product form stationary distribution. We finally present exact tail asymptotics of two 
boundary measures that are associated the stationary distribution. 

Closely related to the first set of our results, Avram et al. [1] studied the variational 
problem that is associated with a multidimensional SRBM. In the two-dimensional case, 
they obtained explicit expressions for the value function I{v) for each point v in the state 
space M^. In particular, they proved that I{v) = mm(I^^\v) , I^'^\v)) , where /^*)(f) is the 
value function of a constrained variational problem with constrained path being restricted 
not to touch boundary 

Fi = {x eRl: Xi = 0} 

of the state space, i = 1,2. To compute /*^*^(f), they showed that an optimal path from 
the origin to v is either a straight line through the interior of the state space or a broken 
path whose first segment travels on boundary F^^i. In either case, they give an explicit 
algebraic expression for I^^\v). Central in their analysis is a notion of "exit velocity" 
that dictates the direction for the broken path to follow when it exits the boundary -F3_j. 
Theorem 3.1 and Lemma 3.2 of this paper give a geometric interpretation of I^^\v). In 
particular, we show that the exit velocity on boundary Fi and the exit velocity on 
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boundary F2 have the following expressions 

respectively, where rf{9), defined in (3.15) in Section 3, is the normal direction of the 
ellipse at point 6 on the ellipse and 6*^*'^'^ is the "symmetry" of 6^"^'^^ on the ellipse; see 
Figure 6 for an illustration of n^{9) and Figure 3 for an illustration of symmetry points. 
Theorem 3.2 gives a geometric interpretation of I{v). Harrison and Hasenbein [ ] re- 
express the solution in [I] for I{v) in an alternative form that is somewhat more explicit, 
and offer further elaboration on the structure of that solution. 

A key step in our analysis is to connect the variational problem with the optimiza- 
tion problems of a linear function over certain convex domains. Those convex domains 
are closely related to the convergence domain of the two-dimensional moment generat- 
ing function of the stationary distribution. Dai and Miyazawa [H] studied the exact tail 
asymptotics of the stationary marginal distribution for the two dimensional SRBM. In 
particular, they characterized the convergence domain through the three geometric ob- 
jects. In this paper, we show that the convergence domain is the intersection of two convex 
domains and the rate function I{v) can be uniquely determined by these two convex do- 
mains. Each of these convex domains can easily be determined by the three geometric 
objects. In this way, we provide a geometric interpretation of the large deviations rate 
function I{v) (see Theorem 3.1 and Theorem 3.2). 

Our second set of results is to characterize geometrically the existence of the product 
form stationary distribution. We prove in Theorem 5.1 that the stationary distribution 
of a two-dimensional SRBM has a product form if and only if 

g[(l,r) ^ g[(2,r)_ 

See Figure 9 for examples. This geometric characterization contrasts with the algebraic 
characterization developed in Harrison and Williams [ 1 ()]. Their characterization, in terms 
of a skew symmetry condition, is valid for an SRBM in an arbitrary dimension. 

Our last set of results (Theorem 6.1) concerns the exact tail asymptotics for two 
boundary measures z/i and z/2 that are associated with the stationary distribution. This 
set of results is closely related to the recent paper [5] in which moment generating functions 
are used to obtain the exact tail asymptotic for stationary marginal distributions. Their 
analysis relies on two asymptotic inversion lemmas to obtain the exact tail asymptotics 
from the singularity of the moment generating functions. Unfortunately, these inversion 
lemmas are not sharp enough for us to study exact tail asymptotics for boundary measures. 
In this paper, we develop a series of three lemmas. Lemmas 6.1, 6.2 and 6.3, to sharpen 
the asymptotic inversion lemmas. These lemmas are likely useful for other problems as 
well. 

This paper consists of six sections. In Section 2, we discuss the variational problem 
of an SRBM and its relation to the large deviations rate function for the stationary dis- 
tribution. In Section 3, we first define three geometric objects and two convex domains 
D*^^) and V'^'^\ We then state the first set of results. Theorem 3.1 and Theorem 3.2. In 
Section 4, we review some results in Avram et al. [•] and prove Theorem 3.1. In Sec- 
tion 5, we give a geometric characterization for the existence of a product form stationary 
distribution. In Section 6, we study exact tail asymptotics of two boundary measures. 
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2 Variational problem and large deviations principle 



In Section 4, we recover results obtained by Avram et al. ['], leveraging recent results 
of Dai and Miyazawa [5]. Those two papers use different sets of notation. For example, 
Avram et al. [1] used F for covariance matrix and 9 for drift vector instead of S and 
in this paper. Furthermore, R in [I] is normalized in such a way that rn = r22 = 1- 
Following [ ], we will not use any normalization in this paper. In Avram et al. [1], the 
inner product (x, y) for vectors x,y ^M."^ is defined as x'^'E~^y where is the transpose of 
X. (All vectors are envisioned as column vectors.) Their inner product certainly simplifies 
computation, but is also potentially confusing because there is no indication of S^^. In 
this paper, we adopt the standard inner product, that is, 

{x, y) = x^y for x, y G 

with the corresponding standard norm ||x|| = a/ (x, x) for each x G M^. 

As in [ ], we assume that S is non-singular and the SRBM data {T,,fi,R) satisfy 

rii > 0, r22 > 0, riir22 - ri2r2i > 0, (2.1) 
^^22/^1 - ri2fJ.2 < 0, rii/i2 - r2ifii < 0. (2.2) 

Under conditions (2.1) and (2.2), the (S, /x, i?)-SRBM is well defined, and its stationary 
distribution exists and is unique; see, for example, [-]. 

We are going to introduce the variational problem associated with the SRBM shortly. 
For that, we first define the Skorohod problem associated with the reflection matrix R. 
Let ©2 = D(]R+,]R2) be the set of functions from M+ to that are right continuous 
on [0, oo) and have left limits in (0, oo). Assume that the reflection matrix R satisfies 
condition (2.1). Then, i? is a completely-iS matrix (for a definition see, for example, the 
second paragraph of the introduction in [']). Therefore, for each x(-) G D^, there exists 
a z{-) G such that z{-), together with some ?/(■) G D^, satisfies the following three 
conditions 

z(t) = x{t) + Ry{t) > 0, t > 0, (2.3) 
y{Q) = and each component of y is non-decreasing, (2-4) 

Zi{t)dy,{t) = 0, 1 = 1,2. (2.5) 



The path z is called an i?-regulation of x. It is not unique in general. Conditions (2.3)- 
(2.5) define the Skorohod problem associated with R. 

Assume conditions (2.1) and (2.2) so that the stationary distribution is uniquely de- 
termined. Let Z{oo) be a random vector whose distribution is the stationary distribution. 
If there is a lower semi-continuous function /(■) from to IR+ that satisfies 

limsup-logP(Z(oo) euB) <- inf /(f), (2.6) 
liminf-logP(Z(oo) e uB) > - inf I{v) (2.7) 

u^oo U vGB° 
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for any measurable B C M^, where B and B° are closure and interior of B, respectively, 
it is said that the large deviations principle (LDP) holds for the stationary distribution, 
and J(f ) is called a large deviations rate function. 

The LDP for the stationary distribution is verified by Majewski [11, Theorem 4] for 
a multidimensional SRBM whose reflection matrix i? is an matrix and the stability 
condition R~^fi < holds. Specializing his result to two dimensions, we see that, under 
condition (2.1), R~^fi < is equivalent to (2.2). Furthermore, again specializing to two- 
dimensional case, for each v G M^, he identifies the rate function I{v) as the value function 
of the following variational problem (VP) 

A^) = inf inf ^ /^(^(t) - /x, S-i(x(t) - /x))dt, (2.8) 

where each a;(-) G H C is an absolutely continuous function on IR+ such that its 
derivative function is locally square integrable with respect to the Lebesgue measure, and 
z{-) G is an i?-regulation of x{-). Dupuis and Ramanan [T] proved an LDP when the 
Ai matrix condition in [11] was relaxed. However, an LDP for a general SRBM remains 
unproven. Despite the lack of such a general LDP, the VP (2.8) is well defined for a 
general SRBM. The focus of this paper is VP (2.8). We use terms value function of the 
VP and large deviations rate function interchangeably even though the latter term does 
not make sense when the corresponding LDP has not been proven. 

For two-dimensional SRBMs, Avram et al. [i] solved VP (2.8) completely. In particu- 
lar, they proved that the value I{v) of (2.8) can be obtained by restricting paths x in one 
of the three types. For each T > and v G M^, define T-L^^^v^T) to be the set of linear 
paths x G from to f such that x{T) = v. Define T-L^^^v, T) to be the set of continuous 
paths X that have two linear segments: during the first segment, an regulated path z of 
X stays on the boundary F2 of the state space and during the second segment, the 
regulated path travels linearly from the boundary to the end point v at time T. Namely, 

n^^\v,T) = !^x en-, x{t) = to for t e[0, s] and x{t) = ^^^{v - z{s)) + so hi 

t G (s, T] for some s G [0, T) and some c G with Ci > and C2 < 

Similarly, define 'H^'^\v,T) to be the set of two segment paths so that during the first 
segment, an regulated path z stays on boundary Fi of the state space. 

Using these sets, let 

/»(t;) = inf , inf ^ T m - fi,J:-\x{t) - ft))dt, z = 0,l,2. (2.9) 

Remark 2.1. We here slightly changed the notation of Avram et al. [1]. Namely, I^^\v) 
and I^'^\v) are P{v) and I^{v) of [1], respectively, while I^^\v) is the same as /''(f) in [ ]. 

We now present the following well known fact (see Theorem 5.1 of Avram et al. [ ] for 
its proof). 
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Lemma 2.1. For an SRBM that satisfies (2.1) and (2.2), the rate function I{v) is obtained 
as 

I{v) = min(/W(f ), I^^\v)), V G M+. (2.10) 
Remark 2.2. I{v) is often written as 

/(t;) =min(/(°)(t;),/«(t;),/(2)(t;)), v eRI- 
However, "H^^) C n^^'^ n "H^^)^ ^nd therefore (2.10) is also vahd. 

Avram et aL [1] exphcitly obtained expressions for I{v) in (2.10) in their Theorem 3.1 
and Lemma 6.1. However, it is difficuh to see how the rate function I{v) depends on the 
SRBM data (S,/i, i?) from their results. In the next section, we use a geometric method 
to compute I{v) directly. In Section 4 we will show that the computation in Section 3 
recovers the results in [1]. 

3 Geometrical view 

Recall that Z{oo) is the random vector whose distribution is the stationary distribution 
of the two-dimensional SRBM. Denote the moment generating function of Z{oc) by if, 
that is, 

^(9) = E (e<^'^(°°)>) , for ^ G (3.1) 

as long as ^p{0) < oo. Dai and Miyazawa [5] identified the convergence domain P, which 
is defined to be 

V = the interior of {9 G M^; ip{9) < oo}. 

In this paper, we provide an alternative characterization of the convergence domain P; 
see Theorem 3.2 below in this section. The characterization allows us to compute the 
value functions I^^\v) and I^'^\v) geometrically for variational problem (2.9) 

3.1 SRBM data and conditions for existence and stability 

To describe the convergence domain V and its alternative characterization, we first intro- 
duce three functions on M^. They are 

^{9) = -^-{9,J:9)-{f^,9), (3.2) 
7,(^) = {R\ 9) = 9^ru + ^2r2i, ^ = 1, 2, (3.3) 

where i?* is again the zth column of the reflection matrix R. These three functions deter- 
mine three geometric objects on the plane M?. Because we assume S is nondegenerate, 
the quadratic equation 7(6*) = determines an ellipse that passes through the origin. 
Two linear equations 71 (^^) = and 72(6') = determine two lines that pass through the 
origin. See Figure 1 for an example. 
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One can check that the transpose of p*^*-* is the ith. row of the matrix det(i?), where 
det(i?) is the determinant of R. Thus, p'-^-' and p*-^-* are orthogonal to the second and first 
columns, respectively, of R. Note that p^^^ = and p*^^^ = p^, where is defined in [1] for 
i = 1,2. Clearly, line 7i(^) = parallels directional vector p^"^^ and line 72(^) = parallels 
directional vector p^^\ It is convenient to assign a direction for each line and from now on 
a directional line is called a ray. We assume ray 7i(^) = points in the direction of p'-^-' 
and ray 72 (^) = points in the direction of p^^\ In Figure 1, these directions are denoted 
by arrows on two lines. 

Specifying the three geometric objects in Figure 1 is equivalent to specifying the SRBM 
data (S, fi,R). A ray passing through the origin can be uniquely determined by its angle 



O2 ii{0) = 62 




Figure 1: Specifying the SRBM data is equivalent to specifying an ellipse and 

two rays, all passing through the origin. The ellipse has normal direction yU at the origin. 
Ray 73_j(6') = has direction p*^*) that intersects with the ellipse at a point 6'^*''"^ 7^ 0, 
i = 1,2. Stability condition (3.4) is satisfied for the left figure and not satisfied for the 
right figure 




Figure 2: The V matrix condition (2.1) is equivalent to that ray 72 (^) = must lie in the 
shaded region 
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in (— 7r, tt], measured counter-clockwise against the 9i axis; the angle is positive if the ray 
points above the axis, and is negative if it points below the axis. Ray one is said to be 
above (below) ray two if angle one is strictly larger (smaller) than angle two. Ray one is 
said to be on ray two if two angles are equal. 

Geometrically, the V matrix condition (2.1) is equivalent to that ray 72 (^) = is below 
ray 71(6*) = 0; see Figure 2 for an example. With the definition of p'-*^ for i = 1,2, the 
stability condition (2.2) can be written as 

(/i,p(i)><0, (/i,p(^)><0. (3.4) 

As discussed in [ ], this is an important observation in the geometrical view of the stability 
condition. Condition (3.4) is equivalent to that the angle between vector /i and each ray 
is more than 7r/2. In Figure 1, the left panel corresponds a stable SRBM and the right 
panel an unstable SRBM. Since fi is the normal direction of the ellipse at the origin, 
condition (3.4) ensures that ray 71(6*) = is not tangent to the eUipse, and therefore it 
must intersect with the ellipse at a point other than the origin. 



3.2 Convergence domain and rate function 

Let 

T = {9e M2;7(^) < 0} 

denote the interior of the ellipse. Its boundary dT is the ellipse itself. Since 7(6') is a 
convex function, F is a convex set. We use 6"-^'™^^^ to denote the right-most point on dT 
and ^(2'™'^^) the highest point; namely, 

Qii,^^^) = aigmaxOi, 1 = 1,2. 

We also define two subsets Fi and F2 of F via 

Ti = {ee F; J3-^{0) < 0} for i = 1,2. (3.5) 

Clearly, Fi is the region inside the ellipse and is below ray 72 (^) = and F2 is the region 
inside the ellipse and is above ray 71 (^) = 0. We use dTi denotes the boundary of Fj, 
z = 1,2. Following [')], we use 9^^'^^ to denote the intersection of ray 72 (^) = with the 
ellipse (other than the origin), and 9^^'"^^ the intersection of ray 71(6') = with the ellipse; 
see Figure 1. We note that 9^^'^'' is proportional to p^^^ and therefore 9^'^^p''^^ = r229^^''^^ . 
Hence, ^{9^^'^^^) = implies that 

(i,r) _ 2r22 (i,r) _ 2ri2 

One can similarly compute 9^''^^ and 9^'^\ Thus, we have 



Condition (2.2) or condition (3.4) is equivalent to 
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92 7iW=0 



^(2, max) 



92 7iW = 




max) 




Figure 3: Three pairs of points: 0^'^'^\ 9^^''^^ and on the elhpse 

02 J)(2,max) ^^^Q^^ ^ Q ^(2,max) 



(l,max) 




( 1 , max) 




Figure 4: An example in which both 9^^'^^ and 9^^'"^^ are lower than 6"-^'™'*^-' 

#''^>0 for 2 = 1,2. 



(3.7) 



Now we define a pair of points 9^^'^^ and 6'*^^'''^ on dT. The point 6'*^^''''' is defined to be 
the "symmetry" of on ^F. When 9^^'"'^ = ^(i--^^^), we define 9'^^''^ to be 9^^'"^^''^ itself. 
When 9^^''^ ^ ^(L^ax)^ ^^^^^ Q(i,r) 

unique point on dV such that 
ef'^)=ef^'^ and ^(^'^V^^^'''^- (3.8) 

Similarly, we can define 9^'^''^\ See Figure 3 for examples of three pairs of points 

and 6'^*'™'*^), i = 1,2. As illustrated in the figure, 9^^'^'' can be higher or lower than 9^^'^\ 



The following is an important condition for the SRBM data: 

^(l,max) ^ 



(3.9) 



In the right panel of Figure 3, condition (3.9) is satisfied, while in the left panel, the 
condition is not satisfied. As it will be explained in Section 4, condition (3.9) is equivalent 
to the fact that the face F2 = {9 & : ^2 = 0} is not reflective, an important term 
introduced in [1]. Under condition (3.9), 9^^'"^"^ is always at most as high as 9^^'^\ Readers 
are warned that condition (3.9) is not equivalent to 



(3.10) 
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The right panel of Figure 4 gives an example that satisfies condition (3.9), but not con- 
dition (3.10). Conditions (3.9) and (3.10) are equivalent if and only if 9^^'"^^ or 9^^'^^ is on 
the right side of ^(2,max)^ namely, 

^(l,r) _ 2(l,r) ^ ^(2,max)_ .3^^^. 



O2 7l{0) = 




O2 7l{0) = 



^(2, max) 



max) 



{e) = o 




max) 



{6)^0 



Figure 5: An example of V^^^ (left) and P*-^^ (right) 

Now we define two subsets V^^^ and P*-^^ of M^. For this, the following points on the 
ellipse dr are convenient: 



^(i,max) 



if ^(^'"1'^=^) ^ dr 

if ^(i,max) ^ QY 



1,2. 



We define 

r'(^) = 1^ G there exists a 0' G F such that 9 < 9' and < ^f'^^}. 



1,2. 



In what follows, we mainly consider T)^^^ because results for V^"^^ are symmetric. Note 
that the definition of V^^^ depends on whether condition (3.9) is satisfied. We have the 
following lemma. 



Lemma 3.1. The domain T)^^^ has the following form 



{9eR^:9< ^(i'"-)} if ^(1---) ^ 9Fi and if 9[''''^ < ^^f '"'""^ 

r^nax n{9eR'':9,< 9\''''>} if ^(i---) ^ dr, and if ^^i'-"^^ > ^^f 



where 



otherwise, 



{6* G M^; there exists a 6*' G F such that 9 < 9'}. 



(3.12) 



(3.13) 



Proof. Assume that condition (3.9) is satisfied. Then, 9^^'^"^ = 9^^'™^^\ and therefore 
1^^^^ = Tmax- We next assume that (3.9) is not satisfied. Then, 9^^''^ higher than 9^^''\ If 
0^^''') is on the left side of ^(2,max)^ ^^1 points 9 G f^^^ are dominated by^6'(^'''); see the left 
panel of Figure 4. Otherwise, 9^^'"^^ must be higher than 9^^'^^^\ and 9^^'"^^ is the lowest 
point on dT fl dV^^^; see the left panel of Figure 5. Thus, (3.12) is obtained. □ 
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Figure 5 illustrates the domains V^^^ and T)^'^\ We are ready to present the first 
theorem on rate function I^^\v). It will proved in Section 4. 

Theorem 3.1. Assume conditions (2.1) and (2.2). For v E M^, 

I^'\v) = sup{{v,e);e eV^'^}, 2 = 1,2. (3.14) 
Corollary 3.1. The value function I^^\v) is a convex function oi v G M^. 

Remark 3.1. Equality (3.14) corresponds to (44) in Theorem 5.1 of for the two- 
dimensional reflecting random walk studied in [2]. 




Figure 6: The normal vector n^{9) that is orthogonal to the tangent line of the ellipse at 

eedv 



For each 6' G M^, define 



(3.15) 



It is easy to see that n^{6) is the outward normal direction of the ellipse a.t 6 G dT. 
Clearly, n^{6) is orthogonal to the tangent line of the ellipse dT at 9. The normal vector 
is illustrated in Figure 6. Since P^*^ is convex and upper bounded, the supremum 

sup{{v,9)-9e 

is finite and is achieved at an extreme point 9 on the boundary of D^*-*. Theorem 3.1 and 
Lemma 3.2 below give a geometric interpretation of I^^\v). 



Lemma 3.2. For each v G M^, 



sup{{v,9);9 e 



(3.16) 



'sup{(t;, 9); 9eT} if ^(i-"^^-) g dTi, 



^sup{(i;, 9); 9eT} if ^ SFi, ^{^'''^ > ^j^''^^"^ is above ^^(^(i''- 



'((2, max) 



The following characterization of the convergence domain D is in a form different from 
the one in [o]. See Figure 7 for an illustration of a convergence domain. Define 



T, = snp{9i>0;9 eVir]V2}, 



1,2. 



(3.17) 
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Figure 7: D = D^^) nD^^) (left), where P^^) and P^^) ^re in Figure 5, and V = P^^) (right) 



Theorem 3.2. Assume conditions (2.1) and (2.2). Then, (a) 

V = v^^'^ nv^'^\ (3.18) 

(b) The large deviations rate function I is obtained as 



in((t;,^(i''^)),(t;,^(2,r)^)^ 



In particular, 



where 



J(e»)=r,, ^ = 1,2, (3.20) 



e(^)=fj) and e^^) = ( J 



Figure 8 gives an illustration the rate function given in (3.19). The following corollary 
is immediate. Inequality (3.21) below is also proved in [■ , Section 8]. 

Corollary 3.2. For each v G 

I{v) >snp{{v,e);e eV}. (3.21) 

Proof of Theorem 3.2. We first prove (a). Here, we adopt the following three categories, 
which are introduced in [5]. 

Category I: '^^ < e\^'^^ and ^^^'^^ < (3.22) 

Category II: ^^^.r) < ^(i,r)^ ^3 23) 

Category III: e^^'^'> < ^^^.H. (3.24) 
It is immediate from the definitions of V^^^ and V^'^^ that 

= G M^; there exists a ^' G F such that 9 < 6', e[ < '^^ and 6'^ < ^^^'^^|. 

12 




Figure 8: The shapes of the rate functions I{v) with \\v\\ = 1 are drawn by thick curves; 
r is inside the elhpse in the right panel and is outside in the left panel 



Hence, 6 e V^'^^ (1 V^^^ implies that 

e < {e^^^^\ef^^Y (3.25) 

for all categories. For Category I, (3.25) and 6 E Fmax implies 9 G V^^^ fl V^'^\ For 
Category II, ^(^.r) < ^(i,r) implies that 6^^'^^ < and therefore ^(2.1^=^^) ^ T2. 

Hence, P^^^ n I^^^) is equal to {6 e R^;9 < 9'-'^''^ by Lemma 3.1. Thus, we have the 
following expression by the symmetry of Categories II and HI. 

f {9 e r^ax; 9 < {9'^'^\ 9^^'^Yy for Category I, 
n p(') = <^ {9eR^;9< for Category II, 

[ {9 eR^;9< ^(^'■')}, for Category HI. 



This implies that 



{9Y\ , for in Category I, 

for Category II, (3.26) 
for Category HI. 



Thus, the r is identical with the one that is obtained in [•')]. Since Theorem 2.1 of [5] says 
that 

^^ = {^Gr„,ax;^<r}, (3.27) 

(3.18) is proved. We next prove part (b). For this, we consider the two cases, r G F and 
r ^ F, separately. First assume that r G F (see, e.g., the left picture of Figure 4). In 
this case, n^(9^^'^^) is above n^(9^'^'^^), and therefore we have the following three cases for 
nonzero v G by Lemma 3.2 and its symmetric version for T)^'^\ If v is below n^{9^'^'^^) 
or on n^(9^'^'^^), then 

j(i)(t;) = (v, P^'^) and I^^\v) = sup{(w, 9);9 E F}. 



13 



If V is between n^^O^"^'^^) and n^(6^^'^^), then 

lW(v) = {v,e^'^'^) and I^^\v) = {v,e^^^'^). 

If V is above n^{9^^'^'^) or on n^(9^^''^^), then 

= sup{(t;, 6); OeV} and I^^\v) = {v, 6^^''^). 

Thus, we have (3.19) for r G F because I^^\v) < I^'^\v) for the first case, and I^'^\v) < 
I^^\v) for the third case. 

We next suppose that r F (see, e.g.. Figure 5). In this case, n^(9^^'^^) is not above 
j7,r^^(2,r)^^ and we similarly have the following three cases. If v is below n^{9^^'^^) or on 
n^{e'^^'''>), then 

lW(v) = {v, 9^'^'^) and I^^\v) = sup{{v, 9)-9e F}. 
If V is between inF [9'^^''^^) and n^(^(^''')), then 

li^)[v) = I^^\v) = sup{(t;, 9)-9e F}. 
If V is above n^{9^'^'^^) or on n'"(^^*^^'''^), then 

/(i)(^) = ^nv{{v,9)-9e F} and /(2)(t;) = {v,9^^''^). 

In the first case, I^'^^v) < I^'^\v), and I^^\v) = sup{^f,0);0 G P}, while in the third 
case I^'^\v) < I^^\v), and I^'^\v) = sup{(f , ^^); G V}. Hence, we always have 

min(j(i)(f ), /(^H^)) = sup{(t;, 0); ^ G I?}. 

Thus, we obtain (3.19) for r ^ F. 

We finally prove (3.20). This is immediate from (3.19), (3.17) and (3.18) if r ^ F. Oth- 
erwise, I{v) = min((w,^(i''")), {v,9^^''^)) by (3.19). This implies (3.20) since {e^^\ 9^^'^^) < 
(e(^), and (e^^), ^(2-'-)^ < ^g^^), ^(i'"-)) for r G F. This completes the proof. □ 



4 Characterization by Avram, Dai and Hasenbein 

Theorem 6.3 of [1] gives explicit expressions for I^^\v) for each v G M^. We will use these 
expressions to prove Theorem 3.1. For concreteness, we will prove Theorem 3.1 for i = 1. 
We first explain how to evaluate 

sup{(t;,e);^ G F} (4.1) 

in the right side of (3.16). 

Recall the normal vector n^{9) for a point 9 G dT. The first component of (^9^^'^^^^) 
is zero and the second component of n^(^(2,max)-j -g ^ero. Namely, the normal at ^(^'™^^) 
is vertical and the normal at 6'^^'™^^) is horizontal. Thus, for any v G M^, there exists a 
unique 9^""'^^ G dT in the segment between 9^^'^^"^^ and 9'^'^'"'^^^ such that n^{9^^'^^) parallels 
V. Clearly, 

{v,9^''^^^) = snp{v,9) . (4.2) 
6»er 

The next lemma gives an analytical expression for {v,9^'"'^^). 
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Lemma 4.1. For each v G M^, 

{v, e^^'^'^) = v/(/i,s-V) {v,^^'v) - {fi, j:-'v) . (4.3) 

Proof. Obviously, the right side of (4.2) is attained by the 6 on the elhpse dT such that v 
parallels the normal n^{6) at the 6. That is, v is proportional to n^{d). Thus, following 
(3.15), we have for some a > 0, 

T,9 + fi = av. 
Substituting 9 = T,^^{av — fj,) in 7(6^) = 0, we have 



a 



Hence, 

{v, 9) = {J:~\av -fi),v) = v/(/i,S-i/i) {v,J:-'v) - (E- V, v) . 
This proves (4.3). □ 

Central in the proofs of Avram et al. ['] is two vectors and defined in (3.2) 
and (3.4), respectively, of their paper. They call the symmetry of around face 
F2 = {z = {zi,Z2) G : Z2 = 0}, namely the nonnegative horizontal axis. (Readers are 
warned again the swap of indexes between [I] and this paper. Anything pertaining to the 
horizontal axis is indexed as 1, with a pair of parentheses, in this paper, and is indexed 
as 2 in [1].) 

It follows from (3.6) and (3.15) of this paper and (3.2) of [1], one immediately has the 
following geometric interpretation 

2 _ „r.z,(i,r) 



a = n 



{9^^'''^). (4.4) 



(Note that our ^(^'""^ in (3.6) has denominator {p^^\ Ep*-^^) because our p^^^ is not normal- 
ized.) To introduce their symmetry using their (3.4) in [1], we introduce 



= ( Jj and = (^Jj . (4.5) 

One can check that the vector (or any other vector) has the following orthogonal 
decomposition under inner product (S^^s,?/) for x,y ^ M^. 

The dominators in (4.6) are due to the lacking of normalization in (4.5). Their (3.4) 
defines via 
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The following equation gives the geometric interpretation of a^. 

^2 = (4.8) 

We leave the proof of (4.8) to the end of this section. 

Following Definition 3.2 and equation (3.5) of [ ], the face F2 is reflective if and only 

if 

al > 0. (4.9) 

The following lemma is immediate. 

Lemma 4.2. Face F2 is reflective if and only if condition (3.9) is not satisfied. 

Now we are ready to prove Theorem 3.1. 

Proof of Theorem 3.1. We prove Theorem 3.1 for i = 1. First assume that condition (3.9) 
is satisfied. In this case, a\<Q and the edge F2 is not reflective. It follows from part (b) 
of Theorem 6.3 in [1], for any v G M^, I^^\v) is given by the expression (3.6) of ['], which 
is the right side of (4.3). Thus, in this case, 

J(l)(^;) = {v.e^"^'^^) = sup{(t;,^^);0 G F} = sup{(t;, ^^); G F^^J, 

which is sup{(i;, 9)]9 E V^^^ because V^^^ = F^ax- 

Now assume that condition (3.9) is not satisfied. In this case, a| > and the edge 
F2 is reflective. If 9'-^''^ is to the left of ^(2,max)^ ^-^^^ ~2 ^ q_ ^^^^^^^ any v G M^, w is 
below a^. It follows from Theorem 6.3 part (a), for any v G M^, I^^\v) is given by 

which is equal to supK^;, 9);9 E V^^^} by Lemma 3.2. Otherwise, we have that 9^^'^^ is to 
the right of or equal to 9^'^'^^^\ Thus, of > 0. It follows from Theorem 6.3, part (b), of 
[1] for any v E that is on or above a^, I^^\v) is again given by the right side of (4.3), 
which is equal to {v,9^'"'^^), and therefore to sup{(f,6');^ G T>^^^} by Lemma 3.2. When 
w G is below a^, it follows again from Theorem 6.3 part (a) that F^'{v) is given by 

{v,9^'^'^), 

which is equal to sup{(f ,6*); 9 E V^^^} by Lemma 3.2. □ 

Proof of equation (4-8). Because 7(^'-^''')) = 0, we have by (4.4) 

{^-\a^ - /i), (a' - /i)> + 2 (/i, S-i(a2 - /i)) = 0, 

and therefore 

(a2,S-ia2> = (/i,/,). (4.10) 

From (4.7), we have 

(a^S-V> = (S-V,e^> = |^^(S-V,e^> = (S-V,e^>. (4.11) 
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Hence, 

^ (S-V,e^)(5^S-V) (a^n^)(a^n^) 
' (S-ie2,e2) (Sn2,n2) 

(S-V,e2)' {a^,n^f 
(S-ie2,e2) ^ (Sn2,n2) 

This and (4.10) imply 

which is equivalent to 7(S~^(a2 — /i)) = 0. Thus, the point 

S-i(a2 -/i) 

is on the ellipse. From equations (4.6) and (4.7), the first component of 

is equal to the first component of S~^(a^ — yu), which is 6^i'^\ Also, from (4.6) and (4.7), 
= a? if and only if {a?,'n?) = 0. The latter is equivalent to that the normal direction 
at is horizontal or 9^^''^ = e^^''^^\ Thus, we have either J:-\a^ - /i) ^ 0^^''^ or 

S"^(a2 - /i) = 6'^^''') = ^(i'™*"). Therefore, we have proved that 

S-i(s2-^) = 0~(i''^), 

which is equivalent to (4.8). □ 



5 Product form stationary distribution 

Harrison and Williams [ I n] proved that a multi-dimensional SRBM has a product form 
stationary distribution if and only if the SRBM data (S, yU, R) satisfies R~^fi < and the 
skew symmetry condition 

2S = i^A^^As + AsA^i/2^ (5.1) 

is satisfied, where, for a matrix A, A^ is the diagonal matrix whose diagonal entries are 
that of A. In this section, we prove the following theorem. 

Theorem 5.1. Assume that (2.1) and (2.2) hold. The two-dimensional SRBM has a 
product form stationary distribution if and only if 

holds, where i?^*'*"^ is the symmetry of 0^^'"^^ defined in (3.8). 
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Figure 9: The locations 9^^''^^ and 9^'^'^^; left panel is for product form; the right panel is 
for non-product form 

Figure 9 gives examples for the locations of 9^^'^^ and 9^"^'^^; left panel illustrates an 
SRBM whose data satisfies (5.2), and the right panel illustrates an SRBM whose data 
does not satisfy (5.2). Equation (5.2) provides a geometric condition on the SRBM data 
for the stationary distribution to have a product form. Of course, the algebraic condition 
(5.1) and geometric condition (5.2) must be equivalent although its verification is not 
obvious. 

Remark 5.1. Avram et al. ['] proved that the skew symmetry condition (5.1) implies 

fi^ = a^ (5.3) 

where is defined in (4.7) and d^ is defined similarly. Because = n^{9^^''^^) and 
d^ = {9^"^'"^^) , where 71^(9) is the outward normal of the ellipse at 9, condition (5.3) 
is equivalent to condition (5.2). Thus, geometric condition (5.2) is necessary for the 
stationary distribution to have a product form. 

We will provide a complete proof of Theorem 5.1 at the end of this section. Before 
that, we provide some background discussion on the product form stationary distribution. 
It is known that the stationary distribution has a density; see, for example, [ ] and [■>]. 
We use Ci^yU) to denote the stationary density of the two-dimensional SRBM. Thus, the 
stationary distribution has a product form if and only if 

C(x, y) = Ci(a;)C2(y) for x,y e M+, (5.4) 

where Q^s are the marginal densities of (. It follows from [10] that when the stationary 
density is of the product form in (5.4), each marginal density d must be exponential. 
Hence, the product form holds if and only if there exist some ai, 0:2 > such that 

C(x, y) = aia2e-("^^+"2^) for x,y>0. (5.5) 

Harrison and Williams [lij] proved that when the stationary density is of the product 
form, a is given by 

a = -2A^^ARR-^fi. (5.6) 
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Associated with the stationary distribution are two boundary measures ui and z/2. The 
measure z/j has support on Fj = {z ^ : Zi = 0}, i = 1,2; see, for example, Section 
2 of [ ] for their definition. It is known that the stationary distribution, together with 
its associated boundary measures, satisfies the basic adjoint relationship (BAR). For a 
statement of BAR in two dimensions, see, for example, (4.1) of [5]. Harrison and Williams 
[ ] proved the necessity of BAR when the reflection matrix i? is an matrix and the 
stability condition R~^n < is satisfied. Assuming the stationary distribution exists, Dai 
and Harrison [3] proved the necessity of BAR for a general reflection matrix R that is 
completely-5. 

Recall the two-dimensional moment generating function f{6) defined in (3.1) for the 
stationary distribution. We denote the moment generating functions of boundary mea- 
sures z/i and z/2 by ipi{02) and ip2{0i), respectively. Then, Dai and Miyazawa [5] derived 
from the basic adjoint relationship the following key relationship among moment gener- 
ating functions: 

'y{eMe) = ^mvM + i2{e)v2{9i) (5.7) 

for any ^ G as long as <^{9) is finite. 

Proof of Theorem 5.1. Remark 5.1 proves the necessity of the theorem. Here, we provide 
a self-contained, alternative proof for the necessity. The technique used in our necessity 
proof will be useful in the sufficiency proof. 

Assume that the stationary distribution has the product form stationary density given 
by (5.5). It follows from [' ] that boundary measure vi has density 

such that 

z/i(M+ X B) = Ci{y)dy for any Borel set B in M+. (5.9) 
Jb 

Similarly, boundary measure z/2 has density 

C2{x) = ^Cix, 0) = ;^aia2e-"^" for x > (5.10) 
2r 22 2r 22 

such that 

U2{A X ]R+) = / C2{y)dy for any Borel set A in IR+. (5.11) 

J A 

It follows from (5.5), (5.8) and (5.10) that 

^(9) = for 9eR^ with 9<a, (5.12) 



ai — ^1 ^2 — O2 

= -^^f^^^, for^2<«2, (5.13) 
2rii(«2 - fc'2) 

¥^2(^1) = ^ ^^^"^"^ ^ for 9, < a,. (5.14) 
2r22(ai - Oi' 
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Plugging these expressions for ip{0), ipi{62), and ip2{0i) into key relationship (5.7), we 
have 

7(^) = ^72(^)(«2 - ^2) + :^7i(^)(«i - ^1) for any 6 < a. (5.15) 
Zr 22 Zr 11 

Since both sides of (5.15) are quadratic functions of 6 and (5.15) holds for infinitely many 
9, it must hold for all 6* e M^. In particular, we have 

7(a) = (5.16) 

That is, a must be on the ellipse dT. Plugging 9 = 9^^'^^ into (5.15) implies that 

«i = ^{1'"-) = ^[i'"^) (5.17) 

because 7(^(^'')) = 0, 72(^^^''')) = and 7i(^(^'')) ^ 0; see Lemma 2.2 of ['] for the latter. 
Plugging 9 = 9^^'^^ into (5.15), we have 



2r 



^72(^~(^'^))(«2-eT^^) = 



(5.1^ 
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because -f{9'^'^'''>) = and ai = 9[^''K If 9'^^''^ = 9^^''\ then 9^^'''> = ^(i'"^^^). This fact, 
together with (5.16) and (5.17), implies that 



a2 — 02 



(5.19) 



If ^(i.'-) ^ 9^^''\ then 72(^(^''')) ^ 0, and (5.18) implies again (5.19). Equations (5.17) and 
(5.19) imply that 

a = 9^^''\ (5.20) 

Similarly, we can prove 

a = (5.21) 
Equations (5.20) and (5.21) imply that (5.2) holds, proving the necessity of the theorem. 



7iW = 



^(2, max) 



(n,-r2) 



^ ( 1 , max) 




2,nax)72W=0 



72(e) = 



71(e) = 



l(2,r)^g(2,r) 




^(l.max) 



Figure 10: The locations of r, 9^^'^^ and 9^'^'^^ for the product form 
We now prove the sufficiency of the theorem. Assume (5.2) holds. Let 



a 



(l,r) _ ^(2,r) 
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^(2, max) 




Figure 11: The case that 9^^'"^^ = ^(i'™^^) for the product form 



Then, a is on the eUipse, and thus 7(0;) = 0. Therefore, 



where 



= /lW(«l-^^l) + /2W(«2-^2), 



(5.22) 



\f2{0)J 2 



S(6' + a) 



Also, 7(0) = gives 



«i/i(0) + a2/2(0) = 0. 
Therefore, there exists a G M such that 



U(o); " 



«2 



(5.23) 



It foUows from (5.22) that 

l{0) = (/i(e)-A;(a2-^2))(«i-^^i)+ (/2(^) + fc(ai-ei))(a2-^^2) 



/l(^)(«l-^l) + /2(^)(«2-^2), 



(5.24) 



where 



/l(^^) = /l(^^)-^(«2-^2), 

Both functions fi{0) and /2(^) are hnear in 9 and (5.23) imphes that /i(0) = and 
/2(0) = 0. If a ^ e^^''\ substituting 6 = O^^'''^ into equation (5.24) yields 



ke^''^) = 

If a = 6'^^'''^ we have 9^^''^'^ = ^(i-™^^), which implies 

d 



(5.25) 



d9. 



-1(0) 



0, 
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from which (5.25) follows again. Thus, we have proved that two linear equations, /2(6') = 
and 72(6*) = 0, both have and 9^^'^^ as their roots. Hence, f2{0) niust be proportional to 
72 (^). Namely, there exists a constant C2 such that 

f^{e) = C2i2{0) for all e e ^5 26) 

Similarly, there exists a constant ci such that 

/i(^) = ci7i(0) for all 6 G (5.27) 
It follows from (5.24), (5.26), and (5.27) that that 7(6') can be written as 

^{9) = ci7i(^^)(ai - 9i) + C272(^^)(«2 - ^2) for all 9 G R\ (5.28) 
Comparing the coefficients of quadratic terms in (5.28), we find 

Q = — , i = l,2. 

Hence, if we define g{9), (71 (6^2) and (72 (6*1) as 

9{0) = n f 2^2 = ^, ^1^1 = w- for ^ < 

«! - 9ia2- 92 a2 -92 ai - 9i 

then g,gi,g2 are the solutions (p,(pi,(p2 to the key relationship (5.7). Let vr be the prob- 
ability measure on with Ci^yV) (5-5) as its density function, and ui and 1^2 be two 
boundary measures defined in (5.9) and (5.11). Thus, we have proved that n, ui and 1^2 
satisfy BAR (4.1) of [ ] for all exponential functions 

with 9 < a. From this, one can argue that BAR is satisfied for all / G C^(]R^), the 
set of functions which together with their ffist and second order partial derivatives are 
continuous and bounded. It follows from [ ] that BAR uniquely determines the stationary 
distribution. Thus, the SRBM must have ({x,y) in (5.5) as its stationary density. □ 



6 Exact asymptotics for boundary measures 

As before, we assume the SRBM satisfies conditions (2.1) and (2.2) so that it has a 
unique stationary distribution. In Section, 5 immediately below (5.6), we introduced two 
boundary measures ui and 1^2 associated with the stationary distribution of the SRBM. 
In this section, we study the exact asymptotics for the tail distribution of i>i. 

The main result of this section is Theorem 6.1, which will be stated shortly. We 
emphasize that Lemmas 6.1, 6.2 and 6.3 that are used in the proof of Theorem 6.1 
constitute a significant contribution of this paper as well. These lemmas can potentially 
to be used for other related problems. For convenience, we write fi{x) ~ f2{x) as a; — )■ cxd 
for two functions /i, /2 on [0, 00) if 

lim fi{x)/f2{x) = 1. 

Recall r defined in (3.17) and the three categories defined through (3.22)-(3.24) in the 
proof of Theorem 3.2. 
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Theorem 6.1. Under conditions (2.1) and (2.2), we have the exact asymptotic: 



hm^^i^^ = 6>0, (6.1) 



for some constant 6 > 0, where k, is given below. 

(a) If Ti < 0^^'^^"^^ ^ then for Category I and Category III, k = 0, and for Category II, 



(b) If Ti = 0^'^^"^^ then for Category I, 



and for Category II, 



^ ~ I _!' 5l(l.r) ^(l,max)' (6-3) 



Qil,r) ^ 5)(l,max) 

^(l,r) ^ ^(l,max)' (6-4) 




We note that, in Category III, one must have ti < 0^^'^^^\ Thus, we do not need to 
consider Category III when ti = 0^^'^^^^\ In Section 5, we introduced the moment gener- 
ating functions (pi{02) and (p2{Gi) for boundary measures z/i and 1^2. Dai and Miyazawa 
[5] studied analytic behaviors of (pi{z) and (p2{z) as functions of complex variable z. In 
particular, they proved that (p2{z) is analytic if ^z < ti and is singular at 2; = ri. Fur- 
thermore, they characterized the nature of the singularities at z = ti. See Lemmas 6.6, 
6.7 and 6.8 of [ ]. We will use the singularity of f2{z) at ti to prove Theorem 6.1. Before 
proving the theorem, we first introduce a series of three lemmas. 

For a nonnegative function / : IR+ — )■ IR+ that is integrable on M_|_, let 

/■oo 

g(z) = / e'''f{x)dx (6.5) 



be its moment generating function. Assume that g is analytic for ^z < a and is singular 
ai z = a for some a > 0. Assuming / is continuous, one can apply the complex inversion 
formula to represent / as a contour integral of g. From this integral representation, one 
then uses the singular behavior of g at z = a to obtain the exact tail asymptotic of /. 
Doetsch [6] is a good reference for the complex inversion formula. Dai and Miyazawa [5] 
have successfully applied this technique to study the exact asymptotic for tail distribution 
of (c, Z(oo)) in an arbitrarily given direction c G M^. In particular, they developed two 
asymptotic inversion lemmas. Lemmas C.l and C.2 of [5], that allow one to directly infer 
the form of exact tail asymptotics of / from the singular behavior of g at z = a. With these 
two lemmas, there is no need to get into contour integral and complex inversion formula. 
For convenience, these two lemmas are reproduced in the appendix of this paper. 

In order to apply one of their asymptotic inversion lemmas, one needs to verify a 
certain set of conditions on the moment generating function g. For example, to apply 
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Lemma B.l, one often needs to verify the sufficient condition (B.2). Namely for some 
positive constants /3, a, b and S such that 



1^(^)1 < for e [0,13] and \^z\ > b. (6.6) 

In the proof of our Theorem 6.1, 

f{x) = z/2(x, oo) for x > 
and the corresponding moment generating function 

g{z) = -{^,{z)~V2m ^z<n. (6.7) 

Whenri < ^J^'"^^"), one can prove that y:>2{z) is bounded in a region given by e [0, ri+e] 
and \'^z\ > b for some e > and some 6 > 0. Still we are unable to prove condition (6.6) 
for some positive 6. Furthermore, we do not know how to check conditions (Bla)-(Blc) 
of Lemma B.l directly. 

To overcome the preceding difficulty, we introduce a new technique. We develop a 
series of three lemmas. These lemmas may have independent interests. Let 



nf)i^)= / fiu)du, x>0, 

J X 

be the tail distribution of /. Inductively, for n > 1, define 

POO 

T-{f){x)= / T^--'\f){u)du, x>0, 



where T^{f) = f ■ Denote the moment generating function of T^{f) by T^{f). Then, it 
is easy to see that 

f{f){z)=z-\g{z)-g{0)), (6.8) 

f\f){z) = z'\g{z)~g{0)-g\0)z), (6.9) 

where g = T°(/) is the moment generating function of / given in (6.5). Thus, it is 
easier to get asymptotics for T"(/) with n > 1 because it is easier for T"^{f){z) to satisfy 
sufficient conditions such as (6.6). Thus, if T"'{f) have the same asymptotics as /, it will 
be more convenient to work with T^{f){z). We will show that / and T'^{f) indeed have 
the same asymptotics if the asymptotic is limited to a certain exponential type. 

We call a function / to be ultimately nonincreasing if there is an xq such that f{x) is 
nonincreasing for x > Xq. We have the following lemma showing the equivalence of tail 
asymptotics. The lemma is a converse of Lemma D.5 of [5]. 

Lemma 6.1. Assume that the nonnegative function / is ultimately non increasing. For 
each integer > 1, if 

(6.1a) T"'(/)(x) ~ bx'^e~"^ as x — )■ oo for some real number k, positive number b and 
nonnegative number a. 
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then f{x) ~ a^hx'^e as x — )■ oo. In particular, (6.1a) implies that T{f){x) ~ 
a"'~^bx^e~°'^ as a: — )• oo for any nonnegative and integrable function / on [0, oo). 

We defer the proof of this lemma to Appendix A. We combine this lemma with the 
classical asymptotic inversion results in Appendix B. The idea is very simple. We just 
use the moment generating function T(/) of (6.8) or T^(/) of (6.9) instead of moment 
generating function g. For this, g of (6.7) is redefined in the following lemma. 

Lemma 6.2. Let 

/■oo 

g{z) = / e'^'uidx) (6.10) 



^0 

be the moment generating function of a finite measure on [0, oo). Assume that g satisfies 
the following two conditions 

(6.2a) there is a complex variable function 

for some positive integers k, some complex number b and positive numbers a, (3 
satisfying < a < /3 such that g{z) — goi^z) is analytic for 'Rz < (3, 

(6.2b) there exists an eo > such that g{z) — go{z) is bounded for all z G C such 
that < ^z < /3 and |2; — a| > cq. 

Then, we have the following exact tail asymptotic for u. 

u(x,oo) ^^—x''~'^e~°"' as X oo. (6.11) 
r(A;) 

The asymptotic of this lemma is less sharp than that of Lemma B.l because the lower 
order term in Lemma B.l has a precise exponential form. However, we do not require for 
i>{x, oo) to be continuous in x. Before proving this lemma, we present similarly a version 
of Lemma B.2. In this case, the asymptotic results are the same, but we need to restrict 
A to be greater than —1. 

Lemma 6.3. Assume u and g as in Lemma 6.2. Assume that there are some a > and 
some 6 E [0, |) such that the following two conditions hold for g 

(6.3a) g{z) is analytic on on Qs{c() = {z E C; z ^ a, | arg(2; — a)| > 6} and is 
bounded on Qs{c() H {z E C; \z — a\ > eo} for some eo > 0, where for 5; G C \ {0}, 
a.Tg{z) G (— vr, vr) is the angle of z; 

(6.3b) for some complex numbers a, b and some real number A satisfying A 7^ and 
A > -1, 

hm ia-z)'{giz)-a)=b. (6.12) 
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Then 



iy{x, oo) ~ as x — )• oo (6.13) 

r(A) 

where, r(A) is the Gamma function as defined in Section 53 of Volume II of [12] (see 
also Theorem 10.14 in its Section 54 for its integral representation), and by convention 
^ = when A = 0, -1, -2, . . .. 

In the reminder of this section, we first prove Theorem 6.1 using Lemmas 6.1, 6.2 and 
6.3. We then prove the latter two lemmas, leaving the proof of Lemma 6.1 to Appendix A. 

The proof of Theorem 6.1 We will use Lemmas 6.2 and 6.3 to prove the theorem. 
Recall that (p2{z) is the moment generating function of boundary measure 1^2 and it will 
play the role of g{z) in these two lemmas. 

(a) For Categories I and III, we have ri = 6^^'^^ < 0^'™^^\ By parts (a) and (c) of 
Lemma 6.6 of [ ], all conditions of Lemma 6.2 with k = 1 are satisfied, and hence the 
lemma implies the theorem with k = for this case. For Category II, if ti 7^ ^1^''^'', then 
ip2{z) has a simple pole at z = ti, and all the conditions of Lemma 6.2 with k = 1 are 
satisfied by parts (b) and (c) of Lemma 6.6 of [ ]. Otherwise, (p2{z) has a double pole at 
z = Ti, all the conditions of Lemma 6.2 with k = 2 are similarly satisfied. In either case. 
Lemma 6.2 implies the theorem with k, being given in (6.2). 

(b) We apply Lemma 6.8 of [ ] because ti = 0^'"""^^^ in this case. We use parts (c) 
and (d) of this lemma for Category I and II, respectively. For Category I, it is not hard 
to see that all the conditions of Lemma 6.3 are satisfied with A = | for 6^'^^ = Q^^''^^^^ 
and A = — I for 6\^'^^ 7^ ^(i.max)^ Therefore the theorem is proved in this case with k being 
given in (6.3). Similarly, for Category II, all conditions of Lemma 6.3 are satisfied with 
A = 1 for d^"''^ = and A = ^ for e'^''^ ^ Thus, the lemma implies the 
theorem with n being given by (6.4). □ 

The proof of Lemma 6.2 Let 

f{x) = u{x, 00) for a; > 0. (6.14) 
We apply Lemma B.l to T{f). Let 

fo{x) = f^^'^'e--^ X > 0, 
b 

9oi^) = (a-z)fe ' zeC\{a}. 

It is clear that the moment generating function of /o is qq. Since go{0) = ba~^ , we have, 
from (6.8), 

T(/o)(.) = z-^ Mz) - ,o(0)) = E 

1=1 ^ ' 
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This and ^o(O) = kha-^^+^^ yield 

f\h){z) = z-^{g,{z)-go{0)-zg'm 
= z-'(T{fo)-g'iO)) 

^1 I a 



V(«- 

k i 



E 



Q,fc+1 ^ — i \^ (Q, _ 2;) 

£=1 



^ ^ ia- z)"" 

e=i m=l ^ ' 

b A,, ^ , a*" 
y \k + 1 - m) 



1 6 ''""^ 



6 

ct"^ (o! — zr Oi'^ (a — -2) 

\ / m=l ^ ' 



We check (Bla), (Bib) and (Blc) of Lemma B.l for f (/) and T'^{h). Note that the 
continuity of / of Lemma B.l is satisfied because T{f){x) is continuous in x for / of 
(6.14). (When Lemma B.l is apphed, T(/) and T'^{fo) serve the roles of g and ^^o; 
respectively, in that lemma.) It is clear that 

fif){z) - f 2(/o)(.) = Z-' {g{z) - g,{z) - {g{0) - ^o(O)) - ^(^'(0) - g'M)) 

has a removal singularity at 2; = 0, and is analytic for ^z < (5 because of condition (6.2a). 
Thus (Bla) is satisfied. By (6.2b) and (6.9), (Bib) is also satisfied. It remains to check 
(Blc). From (B.2) and (6.9), we only need to check that 



00 -1 

iyx 

00 



(3 + ty 



uniformly converges for x > T for some T > 0, but this is already proved on page 237 of 
[6]. Hence, Lemma B.l yields 



Since is a positive integer, this implies 

r(/)(x) ~ A-i-x'^-^e--^ as X ^ 00. 
a'^ 1 [k) 

Applying Lemma 6.1 with n = 1, we conclude (6.11) and thus prove the lemma. □ 

The proof of Lemma 6.3 Let / be defined in (6.14). We first apply Lemma B.2 to 
T(/) to obtain the exact tail asymptotic for T(/). Because of (6.14), we have 

f{f){z) = z-'{giz)-giO)-g'iO)z) (6.15) 
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Obviously, this T{f){z) satisfies the conditions (B2a) and (B2b) by (6.3a) of this lemma. 
Thus, we only need to verify (B2c). 

We consider three cases separately, depending on A > 0, A < or A = 0. We first 
assume that A > 0. In this case, we put = 0, then we have, by (6.12) and (6.15), 

hm {a-z)^fU){z) 

= lim z-\a-z)\g{z)-a + a-g{Q)-g'{Q)z) 
= li_m z~'^(a — z)^[g{z) — aj = ha~'^ . 

We next assume that — 1 < A < 0. Then, we must have, by (6.12), 

a = g{a). 

and therefore it follows from (6.15) that 

T(/)(z) - f (/)(«) = z-'{g{z) - giO) - g'iO)z) - a-'{gia) - giO) - g'{0)a) 

= z-^giz) - a) + ^^^^4%^(^(«) - ^(0)) - —9'iO)- 
^ z^a-^ za 

Hence, letting d = T(/)(a), we have, by (6.12) and the fact that A > —1, 

hm {a-zy{f{f){z)-d)=ba-', 

We finally assume that A = 0. In this case, letting d = a~^(a — g{0) — g'{0)a), (6.12) and 
(6.15) yield 

lim {f{f){z) -d) = a-\b + a - ^(0) - g{0)a) -d = ba'^ 

Thus, (B.3) is verified for T{f) and cq = ba^"^ in all cases. Hence, by Lemma B.2 

T(/)(x) ~ ^^x^-V"^ asx^oo. (6.16) 

(When Lemma B.2 is applied, T(/) serves the role of g in that lemma.) Finally, (6.13) 
and thus the lemma follows from Lemma 6.1 and (6.16). □ 
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A Proof of Lemma 6.1 



We only prove the claim for n = 1 because the case for > 2 is iteratively obtained. 
Since / is ultimately non increasing, we can find Xq > such that f{x) is non increasing 
for X > Xq- By the assumption, we have, for any e > 0, there exists some Xi > Xq such 
that 



(1 - e)bx^e-"'' < T{f){x) = / f{u)du < (1 + e)6a;'^e-"^ Vx > xi. (A.l) 

J X 

Hence, for each 6 > and each x > xi + 6, 

f{u)du < (1 + e)6(x - 5)"e-"(^'^) - (1 - e)6x^e-°^ 



x—S 



X / / \ \ X 

Since f{u) is non increasing in u, we have 



^ - i] te'^e-"" + e(( e"^ + 1 ) bx'^e 



6 lim sup -li^ < (e"^ - 1)6 + e(e"^ + 1)6. 



rr-Kp — ax 



Since e can be arbitrarily small, we obtain 



lim sup < -Je"^ -l)b. 



Letting 5 go to zero, we arrive at 



f{x\ 

lim sup ^ ' < ah. (A.2) 



Similarly, for each 6 > and each x > Xi, 

fX+5 



I f{u)du > (1 - e)6x"e-"" - (1 + e)6(x + (5)"e-"("+'') 



which yields 



1 _ ( ^±l\ " 6x"e-"^ -e[l+ " e""'' ) fex'^e""" 



liminf ^^^^ >ab. (A.3) 



Hence, (A.2) and (A.3) conclude f{x) ~ abx'^e 



B Asymptotic inversion lemmas 

In this appendix, we state Lemmas C.l and C.2 in [ ]. Actually, we here present a refined 
version of Lemma C.l, whereas Lemma C.2 is unchanged. This refinement is closer to 
Theorems 35.1 of Doetsch [(]]. The reason for modifying Lemma C.l is clear in the proof 
of Lemma 6.2 in Section 6. For two functions hi and /12, "we say hi = o(/i2) as x — )■ 00 if 

lim hi{x)/h2{x) = 0. 
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Lemma B.l. Let g be the moment generating function in (6.5) of a nonnegative, con- 
tinuous and integrable function /. Assume the following conditions are satisfied for some 
integer m > 1: 

(Bla) there is a complex variable function 



j 



for some positive integers kj and some positive numbers /3, aj, cf^ for £ = 1, . . . , fc^ 
and j = 1, . . . m with < ai < . . . < am-i < ctm < P such that g{z) — go{z) is 
analytic for "^z < /3, 

(Bib) g{z) uniformly converges to as 2; — )■ oo for < 3^2; < /3, 

(Blc) for some T > 0, e~'^y^g{f3 + iy)dy uniformly converges for x > T. 

Then 

/(^) = E E ^^'"'e-'^^- + o(e-^-), X ^ 00. (B.l) 

Remark B.l. Conditions (Bib) and (Blc) are satisfied if, for some constants a,b,6 > 0, 
\9i^)\ < for ^z e [0,/3] and \^z\ > b. (B.2) 

Lemma B.2. Let / and g be two functions in Lemma B.l. Assume that the following 
two conditions hold for some a > and some 5 G [0, |): 

(B2a) g{z) is analytic on Qs{c() = {z E C; z a, | arg(z — a)\ > 6}, where for 
z G C \ {0}, arg(2;) is the angle of z, 

(B2b) g{z) — )■ as |z| — )■ 00 for z G Qs{c(), 

(B2c) for some complex number d, some A G M and cq G M, 

hm {a-z)^{g{z)-d) =Co. (B.3) 



Then 



/(^) = ^^'-V°^(l + o(l))- (B.4) 
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